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Abstract. We consider nonlinear parabolic evolution equations of the form 
dtu = F (t,x, Du, D^uj , subject to noise of the form H {x, Du) o dB where 
H is linear in Du and odB denotes the Stratonovich differential of a multidi- 
mensional Brownian motion. Motivated by the essentially pathwise results of 
[Lions, P.-L. and Souganidis, P.E.; Fully nonlinear stochastic partial differen- 
tial equations. C. R. Acad. Sci. Paris Ser. I Math. 326 (1998), no. 9] we 
propose the use of rough path analysis [Lyons, T. J.; Differential equations 
driven by rough signals. Rev. Mat. Iberoamericana 14 (1998), no. 2, 215-310] 
in this context. Although the core arguments are entirely deterministic, a con- 
tinuity theorem allows for various probabilistic applications (limit theorems, 
support, large deviations, ...). 



Let US recall some basic ideas of (second order) viscosity theory [6l [7] and rough 
path theory [Ml HZj . As for viscosity theory, consider a real- valued function u = 
u {x) with X G and assume u G is a classical supersolution, 

-G{x,u,Du,D^u) > 0, 

where G is a (continuous) function, degenerate elliptic in the sense that G {x, u,p, A) < 
G {x,u,p, A + B) whenever B > in the sense of symmetric matrices. The idea 
is to consider a (smooth) test function ip which touches u from below at some 
point X. Basic calculus implies that Du (x) = Dip (x) , D^u (x) > Dip (x) and, from 
degenerate ellipticity, 



This suggests to define a viscosity supersolution (at the point i) to G = as a con- 
tinuous function u with the property that (jl.ip holds for any test function which 
touches u from below at x. Similarly, viscosity subsolutions are defined via testfunc- 
tions touching u from above and by reversing inequality in p.ip : viscosity solutions 
are both super- and subsolutions. Observe that this definition covers (completely 
degenerate) first order equations as well as parabolic equations, e.g. by consider- 
ing i9t — F = on R+ X R" where F is degenerate elliptic. The resulting theory 
(existence, uniqueness, stability, ...) is without doubt one of most important recent 
developments in the field of partial differential equations. As a typical result, the 
initial value problem {dt — F) u = 0, u (0, ■) = uq G BUG (R") has a unique solution 
in BUG ([0, T] x M") provided F = F{t, x, Du, D^u) is continuous, degenerate eUip- 
tic and satisfies a (well-known) technical condition (see condition [1] below). In fact. 



1. Introduction 



(1.1) 
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uniqueness follows from a stronger property known as comparison: assume u (resp. 
v) is a supersolution (resp. subsolution) and uq > vq] then u > v on [0,T] x M". 
A key feature of viscosity theory is what workers in the field simply call stability 
properties. For instance, it is relatively straight-forward to study {dt — F) u = 
via a sequence of approximate problems, say {dt — F") u" = 0, provided F"- F 
locally uniformly and some apriori information on the u" (e.g. locally uniform con- 
vergence, or locally uniform boundednestQ- Note the stark contrast to the classical 
theory where ones has to control the actual deriviatives of m". 

The idea of stability is also central to rough path theory. Given a collection 
{Vi,...,Vd) of (sufficiently nice) vector fields on M" and z e ([0,r],M'') one 
considers the (unique) solution y to the ordinary differential equation 

d 

(1.2) y{t)^J2^,{y)7J{t), y(0) = yoeM". 

i=i 

The question is, if the output signal y depends in a stable way on the driving 
signal z. The answer, of course, depends strongly on how to measure distance 
between input signals. If one uses the oo norm, so that the distance between 
driving signals z, 2 is given by 1^ — 5|oo-[o t] ' then the solution will in general not 
depend continuously on the input. 



Example 1. Take n ~ l,d ~ 2, V — (Vi,V2) = (sin (•) , cos (•)) and yo — 0. 
Obviously, 

z" (t) = (-cos (2TTn'^t) , - sin (2Trn^t) 
\n ^ ' n ^ ' 

converges to in oo-norm whereas the solutions to y" — V (y") z", j/q — 0, do not 

converge to zero (the solution to the limiting equation y = 0). 

If |z — zj^.jQ is replaced by the (much) stronger distance 

1^ ~ ^ll-var;[0,T] = ^Up 2^ |^ti,t,+ i " ^ti,t,+ i| , 
{UK10,T]^ 

it is elementary to see that now the solution map is continuous (in fact, locally 
Lipschitz); however, this continuity does not lend itself to push the meaning of 
()1.2p : the closure of (or smooth) paths in variation is precisely VF^'-'^, the set 
of absolutely continuous paths (and thus still far from a typical Brownian path). 
Lyons' theory of rough paths exhibits an entire cascade of (p- variation or 1/p- 
Holder type rough path) metrics, for each p e [l,cx)), on path-space under which 
such ODE solutions are continuous (and even locally Lipschitz) functions of their 
driving signal. For instance, the "rough path" p- variation distance between two 
smooth M''- valued paths z, z is given by 

/ \ ^ I 0) ~{j) 

max sup > \Zf'f.,—Zf'f. 

J=h-M V(t,)c[o,T]^l " " 

where z^^'j = J dz^^ ® ■ ■ ■ ® dzr- with integration over the j-dimensional simplex 
{s < ri < ■ ■ ■ < rj < t}. This allows to extend the very meaning of (|1.2p . in a unique 



^What we have in mind here is the Barles-Perthame method of semi-relaxed limits. We shall 
use this method in the proof of theorem [T] and postpone precise references until then. 
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and continuous fashion, to driving signals which live in the abstract completion of 
smooth M'^-valued paths (with respect to rough path p-variation or a similarly 
defined 1 /p-H61der metric) . The space of socalled p-rough path^ is precisely this 
abstract completion. In fact, this space can be realized as genuine path space, 

where G'^' (W'') is the free step-[p] nilpotent group over M.'^, equipped with Carnot- 
Caratheodory metric; realized as a subset of 1 + tIPl (R'^) where 

is the natural statespace for (up to [p]) iterated integrals of a smooth R'^-valued 
path. For instance, almost every realization of rf-dimensional Brownian motion B 
enhanced with its iterated stochastic integrals in the sense of Stratonovich, i.e. the 
matrix-valued process given by 

(1.3) B'odB^\ 

yields a path B (oj) in (R"^') with finite l/p-H61der (and hence finite p- variation) 
regularity, for any p > 2. (B is known as Brownian rough path.) We remark that 
_ 1^ 1^ B + A where A :— Anti (S*^^^) is known as Levy's stochastic area; in 
other words B (a;) is determined by (B,A), i.e. Brownian motion enhanced with 
Levy's area. 

Turning to the main topic of this paper, we follow [2T| \22\ [23] in considering 
a real-valued function of time and space u = u {t,x) G BUC([0, T] x R") which 
solves the nonlinear partial differential equation 

d 

(1.4) du = F {t,x,Du,D'^u) dt + ^H,{x,Du)dz' 

i=l 

= F {t,x,Du,D^u) dt + H {x,Du)dz 

in viscosity sense. When z : [0, T] — > R"* is then, subject to suitable conditions 
on F and H, this falls in the standard setting of viscosity theory as discussed above. 
This can be pushed further to z G W^'^ (see e.g. [HJ Remark 4] and the references 
given there) but the case when z = z (t) has only "Brownian" regularity (just below 
l/2-H61der, say) falls dramatically outside the scope of the standard theory. The 
reader can find a variety of examples (drawing from fields as diverse as stochastic 
control theory, pathwise stochastic control, interest rate theory, front propagation 
and phase transition in random media, ...) in the articles [551 HH] justifying the 
need of a theory of (non-linear) stochastic partial differential equations (SPDEs) 
in which z in (|1.4[) is taken as a Brownian motiorQ. In the same series of articles 
a satisfactory theory is established for the case of non-linear Hamiltonian with no 
spatial dependence, i.e. H = H {Du). The contribution of this article is to deal 
with non- linear F and H = H (x, Du), linear in Du, although we suspect that the 
marriage of rough path and viscosity methodology will also prove useful in further 
investigations on fully nonlinear (i.e. both F and H) stochastic partial differential 



'In the strict terminology of rough path theory: geometric p-rough paths, 
in which case l|1.4|l is understood in Stratonovich form. 
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equation^ To fix ideas, we give the following example, suggested in [22] and 
carefully worked out in [S] H] . 

Example 2 (Pathwise stochastic control). Consider 

dX = b {X; a)dt + W (X; a)odB + V {X) o dB, 

where b,W,V are (collections of) sufficiently nice vector fields (with b,W depen- 
dent on a suitable control a = a (t) € A, applied at time t) and B,B are multi- 
dimensional (independent) Brownian motions. Define 



V (x, t; B) = inf I 



g{X^/)+ f{X':'\as)ds 



B 



where X^'* denotes the solution process to the above SDE started at X (t) = x. 
Then, at least by a formal computation, 

dv + inf [b (x, a) Dv + LaV + f (x, a)] dt + Dv -V (x) o dB = 

with terminal data v{-,T) = g, and La — in Hormander form. Setting 

u (x, t) = V (x,T — t) turns this into the initial value ( C'auchy) problem, 

du — inf [b (x, a) Du + LaU + / (x, a)] dt + Du ■ V ix) o dBx^- 

with initial data u(-,0) = g; and hence of a form which is covered by theorem\^ 
below. Indeed, H — {Hi,H2), Hi {x,p) — p ■ Vi (x), is linear in p. (Moreover, the 
rough driving signal in theorem]^ is taken as Zt :— B^-t (w) where B (w) is a fixed 
Brownian rough path, run backwards in time^) 



Returning to the general setup of (|1.4p . the results [2Tl[22l[23| are in fact pathwise 
and apply to any continuous path z S C ([0,T] ,M'^), this includes Brownian and 
even rougher sources of noise; however, the assumption was made that H = H {Du) 
is independent of x. The role of x-dependence is an important one (as it arises in 
applications such as example [2]): the results of Lions-Souganidis imply that the 
map 

zeC^ ([0, T] ,R'^)^u (•, •) e C ([0, T] , M") 

depends continuously on z in uniform topology; thereby giving existence/uniqueness 
results to 

d 

du = F {t, X, Du, D'^u) dt + ^Hi (Du) dz' 

i=l 

for every continuous path z : [0, T] — s- R''. When the Hamiltonian depends on x, this 
ceases to be true; indeed, take F = 0,d = 2 and Hi {x,p) = pVi (x) where Vi, V2 
are the vector fields from example [U Solving the characteristic equations shows 
that u is expressed in terms of the (inverse) flow associated to dy = Vi (y) dz^ + 
V2 (y) dz^ , and we have already seen that the solution of this ODE does not depend 
continuously on z = (^z^ , z^) in uniform topologjQ. 



^The use of rough path analysis in the context of nonhnear SPDEs was verbally conjectured 
by P.L. Lions in his 2003 Courant lecture. 

'^Alternatively, the proof of theorem [T] is trivially modified to directly accomodate terminal 
data problems. 

^We shall push this remark much further in theorem [2] below. 
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Of course, this type of problem can be prevented by strengthening the topology: 
the Lyons-theory of rough paths does exhibit an entire cascade of (p-variation or 
1/p- Holder type rough path) metrics (for each p > 1) on path-space under which 
such ODE solutions are continuous functions of their driving signal. This suggests 
to extend the Lions-Souganidis theory from a pathwise to a rough pathwise theory. 
We shall do so for a rich class of fully- nonlinear F and Hamiltonians H {x, Du) 
linear in Du. This last assumption allows for a global change of coordinates which 
mimicks a classical trick in SPDE analysis (Tubaro, Kunita, ...) where a SPDE is 
transformed into a random PDE (i.e. one that can be solved with deterministic 
methods by fixing randomness). In doing so, the interplay between rough path 
and viscosity methods is illustrated in a transpart way and everything boils down 
to combine the stability properties of viscosity solution with those of differential 
equations in the rough path sense. We have the following result^ 

Theorem 1. Let p > 1 and {z'^) C ([0, T] , K*^) be Cauchy in (p-variation) 
rough path topology with rough path limit z e (jo,p-var ^jq^ jij ^ qIp] (jgdj j ^ Assume 

ul e BUG (R") -^uoE BUG (R") , 

locally uniformly as e — > 0. Let F — F {t, x,p, X) be continuous, degenerate elliptic, 
and assume that dt — F satisfies ^^^^^ -invariant comparison (cf. definition[l\ below). 
Assume that V = {Vi, . . . , Vd) is a collection of Lip''^^ (R"; R") vector fields with 
7 > p. Consider (necessarily uniqu^) viscosity solutions G BUC([0, T] x R") 
to 

(1.5) du" = F {t,x,Du'',D^u'') dt- Du" {t,x) -V {x)dz^t) , 

(1.6) u"{0,-) = u"o 

and assume that the resulting family {u" : e > 0) is locally uniformly bounde^. Then 

(i) there exists a unique u G BUG([0,T] x R"), only dependent on z and uq but 
not on the particular approximating sequences, such that u" u locally uniformly. 
We write (formally) 

(1.7) du = F {t,x,Du,D^u) dt- Du{t,x) -V {x)dz{t) , 

(1.8) u{0,-) = uo, 

and also u — when we want to indicate the dependence on z; 

(ii) we have the contraction property 

where u^ is defined as limit of ijP , defined as in 11. 5)] with u" replaced by vP through- 
out; 

(Hi) the solution map (z,uo) ^-^ u^ from 
is continuous. 



'Unless otherwise stated we shall always equip _B(7C-spaces with the topology of locally uniform 
convergence. 

^This follows from the first 5 lines in the proof of this theorem. 

simple sufficient conditions is boundedness oi F {■, ■, 0, 0) on [0, T] X R", and the assumption 
that Uq — » UQ uniformly, as can be seen by comparison. 
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2. Condition for comparison 

We shall always assume that F = F {t,x,p,X) is continuous and degenerate 
elliptic. A sufficient conditior0 for comparison of BUC-solutions to dt = F on 
[0, T] X R" is given by 

Condition 1 ( 6, (3.14)]). There exists a function 9 : [0, oo] — > [0, 00] with 9 (0+) = 
0, such that for each fixed t G [0, T], 

(2.1) F {t, x,a{x - x) ,X) ~ F {t, £,a{x - x) ,Y) < 9 (a\x - xf + \x - 5|) 

for all a > 0, x,x Cz M" and X,Y Cz S" (the space of n x n symmetric matrices) 
satisfy 

-3.(J ?),(- _V).3„(i, -/). 

Furthermore, we require F = F (t, x,p, X) to be uniformly continuous whenever 
p, X remain bounded. 

Although this seems part of the folklore in viscosity theorj0 only the case when 
R" is replaced by a bounded domain is discussed in detail in the literature ([6] (3.14) 
and Section 8] or [Jl Section V.7, V.8]; in this case the very last requirement on 
uniform continuity can be omitted). For this reason and the reader's convenience 
we have included a full proof of parabolic comparison on [0, T] x R" under the above 
condition in the appendix. 

Remark 1 (Stability under sup, inf etc). Using the elementary inequalities, 

|sup (.) — sup (*)! , |inf (.) — inf (*)| < sup |. — *| 

one immediately sees that if F^, Fj^fj satisfy for ^,f3 in some index set with 

a common modulus 9, then inf^ .F!^, sup^ inf^ i^^ ,y etc again satisfy 112. Similar 
remarks apply to the uniform continuity property; provided there exists, for any 
R < 00, a common modulus of continuity an, valid whenever p, X are of norm less 
than R. 

3. Invariant comparison 

To motivate our key assumption on F we need some preliminary remarks on the 
transformation behaviour of 

Du = {diu, . . . , , dnu) , D^M = (diju)^.^^^^^ 

under change of coordinates on R" where u — u {t, ■), for fixed t. Let us allow the 
change of coordinates to depend on t, say v{t,-) := u{t,(f)f (•)) where 0j : R" 



which actually implies degenerate ellipticity, cf. page 18 in [B] (3.14)]. 
^^E.g. in Section 4.4. of Barles' 1997 lecture notes, www.phys.univ- 

tours.fr/~barles/Toulcours.pdf, or section V.9 in [7]. 
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K" is a difFeomorphism. Differentiating v (t, (jj^^ {■)) — u{t,-) twice, foliowed by 
evaiuation at (j)^ (y), we liave, witli summation over repeated indices, 

diU {t, (j>t (x)) = dkv {t,x) di(l)t^'''\^^i^^) 

dijU {t, 0t (x)) = dklV {t,x) 9i(/'t~^'''Ut(:r)5j'^t~^''Ut(a;) + 9kV (t,x) dij(j)t^'''\^^(^^y 

We sfiaii write tliis, somewfiat impreciseij{3 but convenient, as 
(3.1) Du|^^(^) = {Dv\x,D(l)t^\^^(^^-j) , 

D^u\m-.) = (^'i'U,^0r'U,(r.) ®^0t~'U,(r.)) + (^«U,^''/'^'U,(.))■ 
Let us now introduce ^^'^^ as the ciass of all flows of C'^-diffeomorphisms of M", 
(f) ^ {(j)^ : t e [0, T]), such that = Id V0 G ^C") and such that 0^ and (j)'^^ have k 
bounded derivatives, uniformly in t e [0, T]. We say that (n) (f) in iff for 
all multi-indices a with |a| < fc 

da4>{n) da<j>t, da<j>{n)~^ da(j)t^ locally uniformly in [0,T] x M". 

Definition 1 (<i>*^'')-invariant comparison; F"^). Let k > 2 and define F'^ {{t, x,p, X)) 
as 

(3.2) 

F {t, 0t {x) , (p, Dcp^ %^^^)) , {X, D(|)^ ® D(j>^ \(^,)) + {p, D^cp^ 

We say that dt = F satisfies ^'^''^ -invariant comparison if, for every (j) £ <i)('=)j 
comparison holds for BUC solutions of dt — F'^ = 0. More precisely, if u is a 
sub- and v a super- solution to this equation (in viscosity sense, both BUC)and 
u (0, ■)<v (0, •) then u < v on [0, T] x R". 

4. Examples 

Example 3 {F linear). Suppose that a(t,x) : [0,T] x M" ^ M"^"' and b{t,x) : 
[0, T] X R" — > R" are bounded, continuous in t and Lipschitz continuous in x, 

uniformly in t £ [O.T]. If F {t, x,p, X) = Tr a {t,x) a {t,x)'^ X -\- b {t,x) ■ p, then 

^^^^ -invariant comparison holds. Although this is a special case of the following 
example, let us point out that F'^ is of the same form as F with a, b replaced by 



<^'''i*'^)L = cr'„ (i, (x))9^(/)t k = l,...,n;m=l,...,n' 

b't'{t,xf = 6* (i,0t (a:))ai0t~^'''|^^(:,)j +^ (cr'„cr^,ay0t~^'''|^^(y) 



k = l. 



By defining properties of flows of diffeomorphisms, t ^ di(j>^ ^''^\^^(^^), dij4>^ ^''^\c/,_^(y) 
is continuous and the -boundedness assumption inherent in our definition of^^^^ 
ensures that a'^,b'^ are Lipschitz in x, uniformly in t £ [0,r]. It is then easy to 
see (cf. the argument of [7| Lemma 7.1]^ that F'^ satisfiws condition [J\ for every 
(j) g (j)(3) ^ y/jjg implies that <i>'-^-' -invariant comparison holds for BU C solutions of 
dt-F'l' ^ 0. 



^"^Strictly speaking, one should view (^Du, D^uj \ . as second order cotangent vector, the pull- 
back of (^Dv, D^v) \x under (p^^ ■ 
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Example 4 {F quasi-linear). Let 



(4.1) 



F (t,x,p, X) — Tr a {t,x,p) a {t,x,p) X +b{t,x,p) . 



We assume b = b{t,x,p) : [0, T] x R" x R" — > R is continuous, bounded and 
Lipschitz continuous in x and p, uniformly in t £ [0,T]. We also assume that 
a = a {t,x,p) : [0,T] x R" x R" — )■ R"^" is a continuous, bounded map such that 

• (j[t,-,p) is Lipschitz continuous, uniformly in {t,p) G [0, T] x R"; 

• there exists a constant ci > 0, such tha¥^ 

\P~<l\ 



(4.2) 



Vp,g e R" : \cr{t,x,p) -a{t,x,q)\ < Ci- 



l + bl + kl 
for all t e [0,r] and a; G R". 

We show that F"^ satisfies condition{J\for every (j) G <I)('^); this implies that - 
invariant comparison holds for dt — F with F given by To see this we proceed 

as follows. For brevity denote 

p = a{x - x) ,J. ^ D(l)t^\^^i^:',,H. ^ D^(l)t^\^^(^.) 

a. = a{t,(j>t{-),{p,J.)),a.^a.a'^,b.^b{t,(j>t{-),{p,J.)) 

so that 

F'f'it,x,p,X) - TT[a^{{X,J,®J,) + {p,H,))]+b, 
= Tr [ J,a, JJX] + 6, + Tr [a, (p, i?,) ] . 

Hence 

F^ {t, i,p, Y)-F^ {t, x,p, X) = Tr [Jjaj JJF - J.a. Jjx] +65 - 5,+Tr [aj {p, H^) - a, (p, H,)]. 

^ „ ' ' V ' ' 

= :(") 



To estimate (i) note that J^axJ^ — Jx'^x {JxC^xY' ■ The R^" x M^" matrix 

( {Jx^x) {Jx^x)^ JxCTx (JxO'S:)'^ \ 
\ (JxCrx){JxO'x) Jx^xiJxCi) j 



is positive semidefinite and thus we can multiply it to both sides of the inequality 

< 3a 



X 

-Y 



L -I 
-I I 



The resulting inequality is stable under evaluating the trace and so one gets 



Tr 



JxCTx {Jx^x)^ ■ X — JjCTj {Jx^x)^ ■ Y 



< 3aTr 



{Jx^x) {Jx^x^ Jx^x {Jx^x) 



-Jx<^x (JxCTx)'^ + JxO'x {JxCr-x)'^ 



SaTr 



~ 3a \\Jx(^x " Jx<^x\ 



i^Jx^x Jx^x^ {Jx^x Jx^x^ 
2 



(using that Tr [. • ."^] defines an inner product for matrices and gives rise to the 
Frobenius matrix norm \\.\\). Hence, by the triangle inequality and Lipschitzness 



^■^A condition of this type also appears also in [l]. 
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of the Jacobian of the flow (which follows a fortiori from the houndedness of the 
second order derivatives of the flow), 

\\Jx<^x — Jx<y x\\ < \\Jx(yx — Jx<^x\\ + \\Jx<yx — JxCTxW 

< \\Jx\\ \Wx — Cr-xW + ll^^a; — Jx\\ x\\ 

< WJxWlWx - CTxW + C2i<J,(l))\x - X\ 

Since a {t, •, q) is Lipschitz continuous (uniformly in (t, q) £ [0, T] x and 4>f (•) 
is Lipschitz continuous (uniformly in t G [0,T]), we can use our assumption |-^.^[ ) 
on a, to see 

(4.3) ||(Tx — fill < {const) X Ix — i| . 



Indeed, 



= \\a (t, 0j (x) ,p- Jx) - a {t, {x) ,p- Ji)|| 

< \\a {t, 0j (x) ,p- J^) ~ a {t, 04 {x),p- J^)\\ 
+ Ik (t, <l^t {S:),P- Jx) - o {t, (f>t {x),p- Js)|| 

< C2[a,(j))\x-x\+ci- 



l + a\x-S:\{\J^\ + \Js\)' 
and, noting that o (jy^ = Id and sup^j ,j.-)g[Qj.]xRn 1 1 | a; || < C3 implies \\Jx\\ 
\\D(l)i^\^^(x)\\ > I/C3, we have 

a\x - x\\Jx - Jil ^ I _| cia\Jx-Ji:\ 



l + a\x-x\i\Jx\ + \M) - ' ' a\x- S:\i\Jx\ + \M) 

~i C4{a,(l>)\x - x\ 
< \x - '"' 



' \x-S:\ i\J^\ + \Ji\) 
< C5 (cr, 0) |a; - i| . 

Putting things together we have 

\{i)\ < C6 (cr, (p)a\x - . 

As for (ii), we have that, 

\bx~bx\ < \b{t,cj)^{x),{p,J:,))-b{t,(j),{x),{p,J,))\ 

+ \b{t, 0i (f ) , {p, J,))-b {t, 0t (£) , (p, Ji) ) I 

< crib)i\(bt{x)-(btii)\ + \p\\Jx-Jx\) 

where cy (b) is the (uniform in t d [0, T]) Lipschitz bound for b (t, •, •). To get the 
required estimate we again use the regularity of the flow. Finally, for (Hi), 

(Hi) = Tr[ax{p,Hs)-ai{p,Hx)]+Tr[ai;{p,Hx)-ax{p,II^)] 

= Tr [as (p, Hi - H^)] + Tr [(aj - a^) (p, H^)] . 

Using Cauchy-Schwartz (with inner product Tr [. • and p — a {x — x) it is clear 
that boundedness of H and a (i.e. sup^, \IIx\ < 00 uniformly in t [0,T] and sim- 
ilarly for a) and Lipschitz continuity (i.e.\Hx ~ H^l < {const) x |a; — 5:| uniformly 
in t €z [0,r] and similar for a) will suffice to obtain the (desired) estimate 

\{iii)\ < cg X a |a: — x\^ . 

Only Lipschitz continuity of a^ — o'x<j'^ requires a discsussion. But this follows, 
thanks to boundedness of sup^. \(Jx\, from showing Lipschitzness of x i-^ Ux ^ 
a {t,(j)f (x) , (p, Jx)) uniformly in t € [0,T] which was already seen in This 
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shows that F'^ satisfies 12. for any 6 e $(3^ To see that F'^ satisfies condition 
[7] if only remains to see that F'^ {t,x,p,X) is uniformly continuous whenever p, X 
remain bounded. To see this first observe that the flow map (j)^ (x), as function of 
{t,x) e [0,T] X R", is uniformly continuous (but not bounded) while the derivatives 
of the (inverse) flow, given by J., H. above, are bounded uniformly continuous maps 
as functions of t, x. One now easily concludes with the fact the observations that 
(a) the product of BUG function is again BUG and (b) the composition of a BUG 
function with a UG function is again BUG. 

Example 5 (F of Hamilton- Jacobi-Bellman type). From the above examples and 
remark[ll we see that ^^^^ -invariant comparison holds when F is given by 

F (t, x,p, X) — inf ■! Tr (7 {t, x; 7) a {t,x;j)'^ X +6 {t, x; 7) ■ p \ , 

the usual non-linearity in the Hamilton-Jacobi-Bellman equation, and more gener- 
ally 

F {t,x,p,X) = mi ITv (j{t,x,p;-i)(j{t,x,p;-if ■ X + b {t,x,p;T)\ 

7^1 L J J 

whenever the conditions in examples\3[ and^ are satisfied uniformly with respect to 

7 e r. 

Example 6 [F of Isaac type). Similarly, ^^^"^ -invariant comparison holds for 

F {t,x,p,X) — snpinf \ Tt a {t,x; (3,^) a {t,x] P,^)^ X + 6 (i, x; /3, 7) • p > , 
f} 1 I I J J 

(such non-linearities arise in Isaac equation in the theory of differential games), 
and more generally 

a;,p,X) = supinf iXr a {t,x,p; a {t,x,p; P,"ff ■ X +b{t,x,p\ (3,"f)\ 
pill J J 

whenever the conditions in examples\^ and^ are satisfied uniformly with respect to 
(3 € B and 7 S F, where B and T are arbitrary index sets. 

5. Some lemmas 

Lemma 1. Let z : [0, T] be smooth and assume that we are given -bounded 

vector field^ V = {Vi, . . . ,Vd) . Then ODE 

dyt ^ V ivt) dzt, te[0,T] 
has a unique solution flow (of -diffeomorphisms) 4> = 4>^ ^ 

Proof. Standard, e.g. chapter 4 in [13]. □ 

Proposition 1. Let z, V and (p be as in lemma\^ Then u is a viscosity sub- (resp. 
super-) solution ( always assumed BUG ) of 

(5.1) ii {t, x)=F {t, X, Du, D'^u) - Du {t, x) ■ V (x) z (t) 

if and only if v {t, x) :— u {t, (j)^ (x)) is a viscosity sub- (resp. super-) solution of 

(5.2) v{t,x) ^ F'l' {t,x,Dv,D^v) 
where F"^ was defined in 
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In particular, if the vector fields are Lip^, 7 > p + 2, p > 1, then they are also C^-boundod. 
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Proof. Set y — 4>t {x) . When u is a classical sub-solution, it suffices to use the the 
chain-rule and definition of F'^ to see that 

v{t,x) = u{t,y) + Du{t,y) ■ (f>^{x) = u{t,y) + Du{t,y) ■ V (y) it 

< F (t, y, Du (t, y) , D^u {t, y)) = F'>' (t, x, Dv (t, x) , D^v {t, x)) . 

The case when u is a viscosity sub-solution of (j5.1|) is not much harder: suppose 
that {t,x) is a maximum of w — ^, where ^ £ ([0,T] x M") and define i/j e 
C2 ([0, T] X M") by ^ (t, y) = e (i, (y)) . Set y = 0,- (x) so that 

F {t, y, Di; [t, y) , il y)) = F* {i, x, (t, x) , D^^ (i, x)) . 

Obviously, (t, y) is a maximum of u — ip, and since u is a viscosity sub-solution of 
(|5.ip we have 

(i, y) + it, y) V (y) z{t)<F (f, y, i?^ (f, y) , I?V (t, y)) . 

On the other hand, ^ (t, x) ~ ip {t, (f>t (x)) implies <^ {i, x) = ip (t, y)~\-Dip {t, y) V (y) i (i) 
and putting things together we see that 

k {I < (t, X, DC (t, x) , D^C (t, x)) 

which says precisely that w is a viscosity sub-solution of (|5.2p . Replacing maximum 
by minimum and < by > in the preceding argument, we see that if u is a super- 
solution of (|5.ip . then w is a super-solution of (|5.2[) . 

Conversely, the same arguments show that if w is a viscosity sub- (resp. super-) 
solution for ()5.2p . then u{t,y) — v {t,(j}~^ (y)) is a sub- (resp. super-) solution for 
dSIl). □ 

6. Proof of the main result 

Proof. (Theorem [H) Using Lemma [l] we see that (j)^ = <j>^ , the solution flow to 
dy = V (y) dz'^ , is an element of <f> = <i>'^^^. Set F'^ :— F'^' . From Proposition [1] we 
know that m'^ is a solution to 

du= = F {t, y, Du"", D'^u^) dt - Du^ {t, y) ■ V (y) dz" [t) , u" (0, •) = u§ 

if and only if is a solution to dt — F^ ~ 0. By assumption of <i>-invariant 
comparison, 

where are viscosity solution to dt — F"^ = 0. Let (p^ denote the solution flow 

to the rough differential equation 

dy^V (y) dz. 

Thanks to Lip''''''^-regularity of the vector fields (jji^ & ^, and in particular a flow 
of C'^-diffeomorphisms. Set F^ = F'^ . The "universal" limit theorem [26] holds, 
in fact, on the level of flows of diffeomorphisms (see [25] and [TS] Chapter 11] for 
more details) tells us that, since tends to z in rough path sense, 

(j)" (j)^ in $ 

so that, by continuity of F (more precisely: uniform continuity on compacts), we 
easily deduce that 

F^ F^ locally uniformly. 
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From the method of semi-relaxed hmits (Lemma 6.1 and Remarks 6.2, 6.3 and 6.4 
in [6], see also 0) the pointwise (relaxed) limits 

V : = lim sup * w"^ , 

V : = lim inf * -y^ , 

are viscosity (sub resp. super) solutions to dt — = 0, with identical initial data. 
As the latter equation satisfies comparison, one has trivially uniqueness and hence 
V :— V — V the unique (and continuous, since v,v are respectively upper resp. 
lower semi-continuous) solution to 

dtv^F^v, v{0,-) = uoi-). 

Moreover, using a simple Dini-type argument (e.g. [6] p. 35]) one sees that this limit 
must be uniform on compacts. It follows that v is the unique solution to 

(hence does not depend on the approximating sequence to z) and the proof of (i) 
is finished by setting 

{t,x) :^v{t,{cf>^,)-'{x)). 

(ii) The comparison — u^\^.^q t]xR" — l'"o ~ "oIoo R" ^ simple consequence of 
comparison for v,v (solutions to dtv — F^v). At last, to see (iii), we argue in the 
very same way as in (i), starting with 

p2„ _^ pz jQgg^jjy uniformly 

to see that v"" v locally uniformly, i.e. uniformly on compacts. □ 

7. Applications to stochastic partial differential equations 

Applications to SPDEs are path-by-path, i.e. by taking z to be a typical realiza- 
tion of Brownian motion and Levy's area, B (w) = (B, A), also known as enhanced 
Brownian motion or Brownian rough path. The continuity property (iii) of our 
theorem [1] easily allows to identify (|1.7p with z = B (uj) as Stratonovich solution to 
the non-linear SPDE 

du^F {t, X, Du, D^u) dt - Du {t, x) ■V{x)o dB, u (0, •) = wq. 

Indeed, under the stated assumptions the Wong-Zakai approximations, in which the 
Brownian B is replaced by its piecewise linear approximation, based on some mesh 
{O, ^, ^ . . . the approximate solution will convergence (locally uniformly on 
[0, T] X M" and in probabihty) to the solution of 

du^ F (t, X, Du, D^u) dt - Du (t, x) ■ V (x) dB, u (0, •) = uq, 

as constructed in theorem[T] Let us give some applications, typical in the sense that 
they have been studied in great detail in the case of classical stochastic differential 
equations. 

(Approximations) Any approximation result to B in rough path topology 
implies a corresponding (weak or strong) limit theorem for such SPDEs: it suffices 
that an approximation to B converges in rough path topology; as is well known 
(e.g. [m Chapter 13] and the references therein) examples include piecewise linear 
-, mollifier - and Karhunen-Loeve approximations, as well as (weak) Donsker type 
random walk approximations [2]. The point being made, we shall not spell out 
more details here. 
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(Twisted approximations) 

The following results implies en pasant that there is no (classical) pathwise theory 
of SPDEs in presence of spatial dependence in the Hamilonian terms. 

Theorem 2. Let V = {Vi,...,Vd) be a collection of C°° -bounded vector fields 
on M" and B a d- dimensional standard Brownian motion. Then, for every a = 
(ai, . . . , QfAr) G {!,..., d}^, N > 2, there exists (piecewise) smooth approxima- 
tions (z*^) to B, with each z'^ only dependent on {B (t) : t £ I?*^} where (i**^) is a 
sequence of dissections of [0, T] with mesh tending to zero, such that almost surely 

z'^ ^ B uniformly on [0, T] , 

but , solutions to 

du^ = F {t, X, Du^,D'^u^) dt - Du^ {t, x) ■ V (x) dz^ , (0, •) = e BUG (M") , 

( with assumptions on F as formulated in theorem QJj converge almost surely locally 
uniformly to the solution of the "wrong" differential equation 

du = [F (t, X, Du, D^u) - Du {t, x) ■ Va {x)] dt - Du {t, x) ■ V (x) o dB 
where Va is the bracket-vector field given by Va = [Vai, [V'aa, • ■ • [K^jv-u ^on]]] • 

Proof. The rough path regularity of B (oj) implies that higher iterated (Stratonvich) 
integrals are deterministically defined; see [24j First thm.]. Doing this up to level 
N yields a (rough path) Sn (B) and we perturbe it in the highest level, linearly in 
the 

[cai, [ea2, • ■ ■ [eajv„i, eQ„]]] -direction 
of Sn (B) viewed as element in the step-iV free nilpotent Lie algebra. This yields a 
(level- iV) rough path B and we can find approximations (z*^) that converge almost 
surely to B in rough path topology (see [9]). One identifies standard RDEs driven 
by B as RDEs-with-drift (driven along the original vector fields by dB, and along 
Va by dt). The resulting identification obviously holds on the level of RDE flows 
and thus 




The flow idcntiflcation then implies that 

du^F it, X, Du, D'^u) dt - Du [t, x) ■ V {x) dB 

is equivalent to the equation with V (x) dB replaced by V (x) d'B-\-Va {x) dt. □ 

Remark 2. The attentive reader will have noticed that the preceding result also 
holds when the Stratonovich differential odB is replaced by dz for some z £ ([0, T] , K'') ; 
it can then be viewed as result on the effective behaviour of a ( deterministic ) non- 
linear parabolic equations with coefficients that exhibit highly oscillatory behaviour 
in time. 

(Support results) In conjunction with known support properties of B (e.g. [12] 

in p- variation rough path topology or [8 for a conditional statement in Holder rough 
path topology) continuity of the SPDE solution as a function of B immediately 
implies Stroock-Varadhan type support descriptions for such SPDEs. Let us note 
that, to the best of our knowledge, results of this type are new for such non-linear 
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SPDEs. In the linear case, approximations and support of SPDEs have been studied 
in great detail [HI [HI US [ll US] . 

(Large deviation results) Another application of our continuity result is the 
ability to obtain large deviation estimates when B is replaced by eB with e — > 0; 
indeed, given the known large deiviation behaviour of (e_B,e^A) in rough path 
topology (e.g. 19j in p- variation and [11 in Holder rough path topology) it suffices 
to recall that large deviation principles are stable under continuous maps. Again, 
large deviation estimates for non-linear SPDEs in the small noise limit appear to 
be new and may be hard to obtain without rough paths theory. 
(SPDEs with non-Brownian noise) Yet another benefit of our approach is the 
ability to deal with SPDEs with non-Brownian and even non-semimartingale noise. 
For instance, one can take z as (the rough path lift of) fractional Brownian motion 
with Hurst parameter 1/4 < 7? < 1/2 , cf. [5] or [1^, a regime which is "rougher" 
than Brownian and notoriously difficult to handle; or a diffusion with uniformly 
elliptic generator in divergence form with measurable coefficients; see jl2| . Much 
of the above (approximations, support, large deviation) results also extend, as is 
clear from the respective results in the above-cited literature. 



8. Appendix: comparison for parabolic equations 



Let u e BUG ([0, T] x M") be a subsolution to dt-F; that is, dtu-F {t, x, Du, D^u) < 
if u is smooth and with the usual viscosity definition otherwise. Similarly, let 
V e BUG ([0, T] X R") be a supersolution. 

Theorem 3. Assume condition]^ Then comparison holds. That is, 

u (0, •) < « (0, •) =^ u<v on [0, T] x R". 

Proof. We follow [6l Section 8]. Without loss of generality, we may assume that 
dtu ~ F (t, X, Du, D'^uj < —c < and that liuit^T u {t, x) — — oo uniformly in 
X £ i?". We aim to contradict the existence of a point (s,z) G (0,r) x i?" such 
that 

u (s, z) — V {s, z) — S > 0. 
To this end, consider a maximum point (i, x,y) G [0, T) x i?" x i?" of 

(i, x,y)^u {t, x)-v {t, y)-^\x-yf - s (\xf + . 

We first argue that, for small (resp. large) enough values of e and a, the optimizing 
time paramter t G [0, T) cannot be zero. Indeed, assuming i = we can estimate 

(5 — 2e|z|^ — (j){s,z,z) 
< 0(0, £,y) 



sup 

x,y 



uo (x) ~ vq (y) -^\x-yf -s (^\xf + 



passing to the limit e ^ 0, a — s- oo on both sides (using lemma [5] below) then leaves 
us with the contradiction 

< S < sup (mq {x) — Vq {x)) < 0. 
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It follows that t E (0, T). Again, the plan is to arrive at a contradiction (so that we 
have to reject the existence of a point (s, z) £ (0, T) x i?" at which u (s, — f (s, z) > 
0) altogether. To this end, let us rewrite 4> (t, x, y) as 

(j) (t, X, y) = {t, x) - {t, y)- -\x-y\ 

where {t, x) — u (t, x) — e \x\^ and (t, y) = v {t, y) + s \y\^ ■ We may apply the 
(parabolic) theorem of sums [IT, Thm 8.3] at (t, x, y) to learn that there are numbers 
a, b and X,Y & S"- such that 

(8.1) {a,a{x-y) ,X) eP^^+u'' {i,x) , {b, a {x - y) ,Y) e P^'-y" {i, x) 

such that a ~ b — and such that one has the estimate (|2.2p . It is easy to see (cf. 
[SI Remark 2.7]) that (|8.ip is equivalent to 

{a,a{x-y) + 2ex,X + 2el) G P'^'+u{i,x), 

(6, a (i - y) - 2ey, F - 2el) G :P^^~w (t, x) ; 

using that 9tM — F (t, x, Du, D^u) < —c and dtv — F (t, x, Dv, D^v) > (always 
understood in the sense of viscosity sub- resp. super-solutions) we then see that 

a- F (i,x,a{x ~ y) + 2ex,X + 2el) < -c, 

b~ F (£, y,a{x-y) ~ 2ey, Y - 2el) > 0. 

Using a ^ b, this implies 

0<c<F{i,x,a{x-y) + 2ex, X + 2el) - F (£, y,a{x-y)- 2ey, Y - 2el) . 

The last step consists in showing that the right-hand-side converges to zero by first 
sending e — s- and then a ^ oo. (This yields the desired contradiction which ends 
the proof.) If e were absent (e.g. set e = throughout) we would estimate 

F {i, X, a{x-y),X)-F (t, y,a {x - y) ,Y) < (^a\x - yf + \x ~ y|) 

and conclude (using ^6, lemma 3.1]) that 

a \x — y\^ , \x — y\ ^ 'AS a ^ oo 

in conjunction with continuity of 6 at 0+. The present case, e > 0, is essentially 
reduced to the case e = by adding/subtracting 

F{i,x,a{x-y),X)-F (f , y,a{x- y) , Y) . 
It follows that c < (i) -t- {ii) -t- {ii) where 

{i) = \F {i,x,a{x-y) + 2ex,X + 2el) - F {i,x,a{x~y) ,X)\ 
(m) = \F{i,y,a{x-y)-2ey,Y -2el)~F{iy,a{x~y),Y)\ 
{Hi) = 9 {a\x - y\^ + \x - y\^ . 

From lemma [31 we see that {&) p = a[x — y) remains, for fixed a, bounded as 
e ^ 0, (b) 2e \x\ and 2e \y\ tend to zero as e ^ 0, for fixed a, and (c) 

a\x — y^ -\- \x ~ y\ Q as (1/a, e) (0, 0) . 

We also note that (|2.2p implies that any matrix norm of X, Y is bounded by a 
constant times a, independent of e. Combining all this information shows that, 
fixed a, 

lim (i) = lim [ii] = 

£— >0 £— >o 
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while (in) as e ^ and a ^ oo. In summary, 

< c < lim limsup [{i) + (ii) + (ii)] = 

a — *oo ^ ^ 

which is the desired contradiction. The proof is now finished. □ 
Lemma 2. Assume uo,vo e BUG (M"). Then 

Q; 2/2 2"\1 

uo{x) - voiy) - -\x - y\ - e {\x\ +\y\ 1 ^ sup[uo{x) - vo{x)] ass^O,i 



sup 



Proof. Without loss of generality M := sup^, [uq (x) — vo (x)] > 0; for otherwise 
replace uq by uq + 2 |M|. Write M^s for the achieved maximum (at x, y, say) of 
the left-hand-side. Obviously, uq {x) — vq {x) — 2s \xf < M^^s for any x and so 



M < lim inf M„ r. 



(It follows that we can and will consider e (a) small (large) enough so that Ma^e > 
0.) On the other hand, |uo| , \vq\ < R < 00 and so 







< M„,e <2R-^\x-yf -e {\xf + \yf) 



from which we deduce ^\x — yf < 2R, or |f — y| < ^/4R/a. By omitting the 
(positive) penaltity terms, we can also estimate 

Ma,e < Uo {x) - Vo (y) 

= Uo {x) - Vo (x) + avo (y-\/4:R/a^ 
< M + ay, {^/AR/c^ 
where ct^q denotes the modulus of continuity of vo- It follows that 

lim sup Ma,e < M 



a — ^00 



which shows that the limMa^g (as e ^ 0, a — > 00) exists and is equal to M. □ 

Lemma 3. Consider a maximum point (i, x, y) € [0, T) x R" x R" of 

<j) {t, x,y)=u {t, x)-v {t, y)-^\x-yf -s (\xf + \yf^ 

where u,v G BUG ([0, T] x M"). Then 

lim sup a{x — y) = < 00, 

lime(|f| + |y|) = 

€—*0 

-\x-y\ +\x-y\j = 0. 

a— »oo 

Proof Without loss of generality M := supjQ yjxi?" (^i^;) — v {t,x)] > 0; for oth- 
erwise replace u by u + 2 |M|. Write Mq, ^ for the (achieved) maximum at some 
i, X, y. Obviously, u {t, x) — v (t, x) — 2e \x\'^ < Ma^e for any {t, x) G [0, T] x i?" and 
so 

M < hm inf M„ 5. 

£— >0 ' 
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(It follows that we can and will consider e (a) small (large) enough so that Ma,e > 
0.) On the other hand, supjQ -pj^Rn \u {t,x) \ , supp ^R" 1" (^i^;)] < R < oo and so 

< < 2i? - I |x - y|' - £ + lyf) . 

(This allows already to deduce 

^\x-y\\ e{\x\' + \y\') <2R 

which readily implies a |a; — y| < 4ai? and e , e |y| < \/2Re.) Let us proceed 
as in the proof of lemma [5] and estimate 

Ma,e < u (£, x)~v {i, y) 

~ u (i, x^ — V (i, x) + (7v (^\/ 4i?/a^ 

< M + ay {^y/m/aj 
from which we see that Ma.e ^ M as e — > 0, a ^ oo. Necessarily then, 
u {i, x)~v {i y) - M^,, = ^\x~y\^ + e {\x\^ + ^ 
as e — > 0, a — > oo. In particular, 

lim ^ \x-y? -> 0. 

E— »o z 

a. — *oo 

Since |x — i/| ^ as e ^ 0, a ^ oo is clear from |a; — y| < AR/a the proof of the 
last estimate is also finished and the proof of the lemma is complete. □ 
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